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Preface 



This thesis for the University Degree of Master of Science is an account of my research at the Insti- 
tute of Physics at the; Autonomous University of San Luis Potosi, Mexico, during the time period 
September 2005 — December 2006. It is devoted to the exact solutions of the Schrodinger equation 
with the hyperbolic Scarf potential (Scarf II) which finds various applications in physics ranging 
from soliton physics to electrodynamics with non-central potentials. In the current literature the 
exact solutions to Scarf II are written in terms of Jacobi polynomials of purely imaginary argu- 
ments and parameters that are complex conjugate to each other. The fact is that the above Jacobi 
polynomials arc proportional to real orthogonal polynomials by the purely imaginary phase factor, 
(— i)", much like the phase relationship between the hyperbolic-and the trigonometric functions, 
i.e. sin(za::) = zsinh(x). These real polynomials, to be referred to as the Romanovski polyno- 
mials, have been largely ignored by the standard mathematical textbooks on polynomials and 
mathematical methods of physics texts although they are required in exact solutions of several 
physics problems ranging from supersymmetric quantum mechanics and quark physics to random 
matrix theory. It is one of the virtues of the present thesis to draw attention to the Romanovski 
polynomials and their importance in the physics with Scarf II. 

I begin with reviewing the five possible real polynomial solutions of the generalized hypergeo- 
nietric differential equation. Three of them are the well known classical orthogonal polynomials of 
Hermite, Laguerre and Jacobi, but the other two are different with respect to their orthogonality 
properties. The family of polynomials which exhibits finite orthogonality (meaning that only a 
finite number of them are orthogonal) are the Romanovski polynomials. 

Next, I solve the one-dimensional Schrodinger equation with Scarf II in terms of the Ro- 
manovski polynomials. Then I focus on the problem of an electron within a particular potential 
that appears non-central in the polar angle coordinate while preserving the rotational invariance 
with respect to the azimuthal angle. I report the new observation that the (one-dimensional) 
Schrodinger equation with the hyperbolic Scarf potential defines the polar angle part of the re- 
spective wave functions. The latter define new non-spherical angular functions. I furthermore 
establish a new non-linear relationship between the Romanovski polynomials to the associated 
Legendre functions and employ it to obtain a non-standard orthogonality integral between infinite 
scries of polynomials. This infinite orthogonality does not contradict the finite one because in the 
latter case the parameters change with the polynomial degree and same does the associated weight 
function. This circumstance allows to satisfy the orthogonality condition for an infinite number 
of polynomials. Finally, I also solve the Klein-Gordon equation with scalar and vector potentials 
of equal magnitudes and given by Scarf II. I conclude that the Romanovski polynomials are the 
most adequate degrees of freedom in the mathematics of the Scarf II potential. 
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Chapter 1 

Introduction. 



The exactly solvable Schrodinger equations (SE) occupy a pole position in quantum mechanics 
in so far as most of them directly apply to relevant physical systems. Prominent examples are 
the quantum Kepler-, or, Coulomb-potential problem and its importance in the description of the 
discrete spectrum of the hydrogen atom [T], the harmonic-oscillator-, the Hulthen-, and the Morse 
potentials with their relevance to vibrational spectra [2] , [3] . Another good example is given by the 
Poschl- Teller potential [4] which appears as an effective mean field in many-body systems with S- 
interactions [5] . Knowing the exact SE solutions is furthermore of interest in testing approximative 
iteration methods [5]. 

There are various methods of finding the exact solutions of a Schrodinger equation for the 
bound states, an issue on which we shall focus in the present work. 

• The traditional method [7| consists in reducing SE by an appropriate substitution of the 
variables to that very form of the generalized hypergeometric equation [8] whose solutions 
are polynomials, the majority of which, and especially the classical ones, being well studied. 

• The second method suggests to first unveil the dynamical symmetry of the potential problem 
and then employ the relevant group algebra in order to construct the solutions as the group 
representation spaces 

• Finally, there is also the most recent and powerful method of the super-symmetric quantum 
mechanics (SUSYQM) which considers the special class of Schrodinger equations (in units 
of ?i = 1 = 2m) that allows for a factorization according to [12] , [12] 

(iJ(z) - e„) Vn(^^) = (^-^ + v{z) - Bn^ iJn{z) = , 

H{z) = A+iz)A-{z) + eo, 

= (±± + Uiz)^ . (1.1) 

Here, H{z) stands for the (one-dimensional) Hamiltonian, U{z) is the so called super- 
potential, and the operators (z) ladder between neighboring solutions. The super-potential 
allows to recover the ground state wave function, ipgst{z), as 

i,g,t{z)-e'i'^^y^''y. (1.2) 

The excited states are then built up on top of 'tpgstiz) through the repeated action of the 
A+(z) operators. The supersymmetric quantum mechanics manages the family of exactly 
solvable potentials presented in Table 1.1. Though this Table contains 10 potentials, the 
variety is in fact not as big. The reason is that these potentials split into two classes 
such that the potentials belonging to same class can be mapped onto each other by means 
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of appropriate point transformations of the coordinates (so called "shape invariance" , see 
Ref. [14l|T6] and references therein). The Coulomb, oscillator, and Morse potentials belong 
to the first class, while the second class contains the Rosen-Morse, Eckart, Poschl- Teller and 
Scarf potentials. There are more exactly solvable potentials. The very popular in nuclear 
physics Woods-Saxon potential, 

r-Rp 

where Vq, Roland a are constant parameters, has been resolved recently in Refs. [15] in terms 
of the Jacobi polynomials. 

• The above methods refer to potentials with unconstrained parameters. However, there are 
also potentials which are exactly solvable only when the parameters obey certain constraints. 
Take as example the sextic anharmonic oscillator, ax^ + bx'^ + cx^ with a > which allows for 
exact solvability [T7] only if the parameters are restricted by the condition, j^ + c~ 0. 

In this case, any wave function corresponding to a level of energy E acts as the generating 
function of a set of polynomials in the energy variable which differ from the solutions of the 
generalized hypergeometric equation. For the first time such polynomials have been reported 
by Bender and Dunne [18j while solving the one dimensional {ID) Schrodinger equation with 
the potential, x^ — (4s — 4 J — 2)a:^, for J a positive integer and s = 1/4, or s = 3/4. The 
two choices of s refer to even- and odd-parity of the eigen-wave functions, respectively. One 
can find a "weight function" with respect to which the Bender-Dunner polynomials appear 
"orthogonal" though the "weight function" neither needs to be positive. 

While the non-rclativistic Schrodinger equation is exactly solvable for the potentials listed 
in Table 1.1 (and for a few more), exact solutions to the relativistic Klein-Gordon- and Dirac 
equations are quite scarce. Besides the relativistic Coulomb problem, the Dirac equation seems 
so far to be amenable to an exact solution only for the relativistic oscillator, a result due to 
Moshinsky and Szczepaniak [19j . Klein-Gordon equations with scalar and vector potentials of 
equal magnitudes reduce to Schrodinger equations and allow for exact solvability (see Chpt. IV 
below). Finally, very recently, low-power potentials of the type, ~ with —l^fi^ —2, 

have been shown to be quasi-solvable in the sense that only the zero energy solutions are exact 
|20j . The present study is devoted to the first item above. 



1.1 The goal. 

Table 1.1 shows that the hyperbolic-Scarf- and the trigonometric Rosen-Morse poten- 
tials (termed to as Scarf II, and Rosen-Morse I, respectively) are solved in terms of 
Jacobi polynomials of imaginary arguments and parameters that are complex conjugate 
to another. 

The goal of this thesis is to solve the Schrodinger equation with the hyperbolic Scarf 
potential [Tl], [lT]-[53] anew and to make the case that it reduces in a straightforward 
manner to a particular form of the generalized real hypergeometric equation whose 
solutions are given by a finite set of real orthogonal polynomials. In this manner, 
the finite number of bound states within the hyperbolic Scarf potential is brought in 
correspondence with a finite system of orthogonal polynomials of a new class. In due 
course, various new properties of the above polynomials are encountered. 

These polynomials have been discovered in 1884 by the English mathematician Sir Edward John 
Routh and rediscovered 45 years later by the Russian mathematician Vsevolod Ivanovich 
Romanovski in 1929 25j within the context of probability distributions. Though they have been 
studied on few occasions in the current mathematical literature where they are termed to as 
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Potential 


Variable 


Wave Function 


Shifted oscillator: ^oj^ {x - |^)' - | 




e^p{-ky^) Hn{y) 


3D-oscillator: ^w^r^ + (l + I) ^ 


1 2 


y^^O+i)exp(-iy) L^^^y) 


Coulomb: + ^) + 


y ^ n+l+l 


y'+iexp(-iy)4-]^l(y) 


Morse: + B'^ exp(-2Q!a:) 
-2B {A + f ) exp(-ax) 


y exp{-ax) 


y-nexp (-iy) L^^^^-"»(y) 


Scarf II : + [B'^ - - ^a)secli^(ax) 
(Poschl- Teller 1) +B{2A + Q;)secli(aa::)tanh(Q;a::) 


y = sinh(ax) 


z"(l + y^y^ exp(-Atan"^(y)) 
xp!l''—'--''—'^\^y) 


Rosen-Morse 11 : A^ + ^ 
-A{A + a) sech2(ax) + 2Btanh{ax) 


y = tanh(ax) 


(l-y)*(l + y)*Pi^-^^^(y) 


Eckart: A"^ + ^ - 2Bcotli(ax) 
+A{A — a)cscli^(Q;x) 


y = coth(ax) 


(y-l)*(y+l)*Pi^-^^^(y) 


Scarf I: -A + [A^ + {A + B'^ ~ Aa) sec^{ax) 
+A{2A — a) tan(aa:) sec(ax) 


y = sin(aa;) 


(l-y)^(l + y)^ 


Poschl-Teller 2: A^ + (B^ + A^ + Aa) cscli2(Q!a;) 

—B{2A + a) coth(aa:)csch(Q!a;) 


y — cosh(Q;a;) 


(y-l)V(y + l)-^ 
xpi'—^^''-'-^\y) 


Rosen-Morse I: A{A — a) csc^(aa::) 
+2Bcot(ax) -A^ + ^ 


y = cot(ax) 


j_n^y2 _|_ 1)-^ exp(acot~-^(y)) 
xPi"^~"-"'-'-"+")(iy) 



Table 1.1: Exactly solvable shape invariant potentials in one dimension (according to [131 [H]). 
All potential parameters are supposed to be positive. Other notations are: s = A/ a, A = B/a, 
a = A/(s — n), si — s — n + a, S2 ~ s — n ~ a, S3 = a — n — s, and S4 — — (s + n + a). The 
orthogonal polynomials of Hermite, Laguerre, and Jacobi have been denoted in the standard way 
as H{x), lI^\x), and Pn^'^\x), respectively. 
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"finite Romanovski" [Hj-El], or, "Romanovski-Pseudo-Jacobi" polynomials [501 j they have been 
completely ignored by the textbooks on mathematical methods in physics, and surprisingly enough, 
by the standard mathematics textbooks as well [5], [21]- [M|- The notion "finite" refers to the 
observation that for any given set of parameters (i.e. in any potential), only a finite number of 
polynomials (finite number of bound states) appear orthogonal. 

The Romanovski polynomials happen to be equal (up to a phase factor) to that very Jacobi 
polynomials with imaginary arguments and parameters that are complex conjugate to each other 
[35|-|37j. much like the sinh(z) = zsin(zz) relationship. Although one may ( but has not to ) deduce 
the local characteristics of the latter such as generating function and recurrence relations from 
those of the former, the finite orthogonality theorem is qualitatively new. It does not copy none 
of the properties of the Jacobi polynomials but requires an independent proof. For this reason, in 
random matrix theory |39j the problem on the gap probabilities in the spectrum of the complex 
circular Jacobi ensemble is preferably treated in terms of the real Cauchy random ensemble, a 
venue that conducts once again to the Romanovski polynomials because, as it will be shown 
below, the weight function of the Romanovski polynomials for certain values of the parameters 
equals the Cauchy distribution. 

The thesis contributes two new examples to the circle of the typical quantum mechanical 
problems [40j . The techniques used here extend the teachings on the Sturm-Liouville theory of 
ordinary differential equations beyond their textbook presentation. 

Before leaving the introduction, a comment is in place on the importance of the Scarf II poten- 
tial. It finds various applications in physics ranging from electrodynamics and solid state physics 
to particle theory. In solid state physics Scarf II is used in the construction of more realistic 
periodic potentials in crystals 0Tj than those built from the trigonometric Scarf potential (Scarf 
I) [42]. In electrodynamics Scarf II appears in a class of problems with non-central potentials 
[331 133] (also see section IV below for more details) . In particle physics Scarf II finds application 
in studies of the non-perturbative sector of gauge theories by means of toy models such as the 
scalar field theory in (1-1-1) space-time dimensions. Here, one encounters the so called "kink -like" 
solutions which are no more but the static solitons. The spatial derivative of the kink-like solution 
is viewed as the ground state wave function of an appropriately constructed Schrodinger equation 
which is then employed in the calculation of the quantum corrections to first order. In Ref. |45j it 
was shown that specifically Scarf II is amenable to a stable renormalizable scalar field theory. 

The thesis is organized as follows. 

• In the next chapter I first highlight in brief the basics of the generalized hypergeometric 
equation, review the classification of its polynomial solutions in various schemes, and con- 
struct for completeness of the presentation, the polynomials of Hermite, Laguerre, Jacobi, 
Romanovski, and Bessel. 

• In chapter III I bring examples for potentials whose exact solutions are given in terms of the 
polynomials presented in the previous chapter. 

• The original contribution of thesis is presented in chapter IV. There I 

— first provide the exact solution of the lZ3-Schrodinger equation with Scarf II in terms 
of the finite Romanovski polynomials. 

— Next I consider the problem of an electron within a particular non-central potential 
along the line of Ref. [321 133] • I made the new observation that the polar angle part 
of the electron wave function in this problem happens to be defined by Romanovski 
polynomials whose parameters depend on the polynomial degree. Within this context, 
the Romanovski polynomials act as designers of new non-spherical angular functions. 

— I explicitly construct and display graphically the lowest five non-spherical-functions 
and compare them with the standard Y{^{9, ip) harmonics. 
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— I furthermore establish a new non-Hnear relationship between the Romanovski poly- 
nomials to the associated Legendre functions and employ it to obtain a non-standard 
orthogonality integral between infinite series of polynomials. This infinite orthogonality 
does not contradict the finite one because in the latter case the parameters change with 
the polynomial degree and same does the associated weight function. This circumstance 
allows to satisfy the orthogonality condition for an infinite number of polynomials. 

— Before closing by a brief summary, I consider a Klein-Gordon equation with scalar and 
vector potentials of equal magnitudes which are given by Scarf II. 
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Chapter 2 

Generalized hypergeometric 
equation. Polynomial solutions. 



All classical orthogonal polynomials appear as solutions of the so called generalized hypergeometric 
equation (the presentation in this section closely follows Ref. [ISj), 



'^{x)y'^{x) + T{x)y[^{x) - A„?/„(x) = , 
ct{x) = ax^ -\-bx + c, ^{x) = xd + e , A„ = n{n — l)a + nd . 



There are various methods for finding the polynomial solution, here denoted by 

Unix) = Pn 



d e 
a b c 



The symbol P„ 



d e 



(2.1) 
(2.2) 



(2.3) 



X makes the equation parameters explicit and stands for a polyno- 



a b c 

mial of degree n, A„ being the eigenvalue parameter, and n — 0,1, 2, 



2.1 Classification scheme of Koepf-Masjed-Jamei. 

In Ref. [2H] the solutions to Eq. (|2.ip have been classified according to the five parameters a, b, c, 
d, and e. Furthermore, a master formula for a generic monic polynomial solutions, P„, has been 
derived by Koepf and Masjed-Jamei, according to them one finds 



P. 



d e 
a b 



E 

k=0 



G["\a,b,c, d, e)x'' 




2ae.—bd 
2aV b'-'-iac 



2a 



2Vb^ - Aac 



b + VlP 



Aac J 
(2.4) 



The a = case is handled as the a ^ limit of Eq. (|2.4p and leads to the appearance of 2-F0 
in place of 2F1 (see Ref. [53] for details). Though the original derivation of this result is a bit 
cumbersome, its verification with the help of the symbolic software Maple is straightforward. 
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2.2 Classification according to Nikiforov-Uvarov. 



In the method of Nikiforov and Uvarov 8 the (polynomial) solutions to the hypergeometric 
equation (12. ip are classified according to the so called weight function, w{x), and built up from 
the Rodrigues formula, 

<^{x)yn{x) + yi{x)y'„{x) - XnVnix) = 0, (2.5) 
2/1(2^) = ;^^^[w(x)cr(x)], Xn = -n{y[{x) + -{n - l)a" (x)) , (2.6) 

with A^i being the yi{x) normalization constant. The weight function is calculated from integrating 
the so called Pearson differential equation, 

^ {a{x)w{x)) = yi{x)'w{x) , (2.7) 

and then plugged into the Rodrigues formula to generate the polynomial solutions as 

N rl^ 

y„(x) = ^— K^)a"(x)], (2.8) 
w{x) ax"- 

with Nn being a normalization constant. The 2/„(x)'s are normalized polynomials and are orthog- 
onal with respect to the weight function w{x) within a given interval, [/i,?2], provided 17(0;) > 0, 
and z/;(x) > holds true for all x S [^1,^2], and a{x)'w{x)x''\x=ii = cr{x)'w{x)x''\x=i2 = for any / 
integer, and with suitable normalization constants, 

h 

w{x)yn{x)yn' ix)dx ^ Snn' , V n, n' G {0, 1 , 2, . . .} . (2.9) 

L 

The Nikforov-Uvarov method combines well with the Koepf-Masjed-Jamei scheme. Indeed, the 
solution to the Pearson equation can be cast into the form 



LU{X) = W 



, , {d~2a)x+{e-b) . 
( / ax^ + bx + c 1 • ^'-''^ 



It shows how one can calculate any weight function associated with any parameter set of interest 
(the symbol used for the weight function makes again the equation parameters explicit). The best 
strategy in recovering the weight function corresponding to a given set of parameters is to consider 
its logarithmic derivative, ^^^j^, and then match the parameters accordingly. In order to illustrate 
this procedure, we pick up one of the examples given in |28j and consider the weight function 

Lo{x) = i-x^ + 3x~ 2)1" = (1 - x)^"{2 - x)i° . (2.11) 

The calculation of the logarithmic derivative gives 

uj'ix) -20x + 30 _ id~2a)x+{e-b) 



uj{x) —x"^ + 3x — 2 ax^ + bx 



(2.12) 



which matches with the parameters a — —1, b — 3, c — —2, d ~ —22, and e = 33. Therefore, the 
resulting weight function is 



uj{x) = W 



-22 33 
-1 3 -2 



(2.13) 



11 



In the notations of Koepf-Masjed-Jamei the Rodrigues formula looks like 



Pn 



d e 
a b c 



= nlZiid+in + k-2)a)Pn 



(ax^ +bx + c)">V 




d e 
a b c 





' d 


e 








■) 




a 


b c 





(2.14) 



The polynomials associated with the weight function in Eq. (|2.1ip can be constructed explicitly 
by Eq. (j2.14p and treated as independent entities and without even knowing that they are no more 
but Jacobi shifted to the interval x Cz [1,2]. The great appeal of combining the master formulas in 
the respective Eqs. (|2.4p . and (|2.14p is that they allow for the direct and pragmatic construction 
of all the polynomial solutions to the generalized hypergeometric equation. 



One identifies as special cases 

• the Jacobi polynomials with a = —1, 6 = 0, c = 1, = —7 — (5 — 2, and e = —7 + S, 

• the Laguerre polynomials with a — 0, b = I, c = 0, d = —I, and e — a + 1, 

• the Hermite polynomials with a = b = 0, c=l,d = —2, and e = 0, 

• the Romanovski polynomials with a = I, b = 0, c = 1, d — 2(1 — p), and e = q with p > 0, 

• the Bessel polynomials with a ^ 1, b = 0, c — 0, d = a + 2, and e — /3. 

These parametrizations will be referred to as " canonical" . Any other parametrization can be 
reduced to one of the above sets upon an appropriate shift of the argument. 

The first three polynomials are the only ones that are traditionally presented in the standard 
textbooks on mathematical methods in physics such like [3T]-[31], while the fourth and fifth seem 
to have escaped due attention. Notice, the Legendre-, Gegenbauer-, and Chebychev polynomials 
appear as particular cases of the Jacobi polynomials. The Bessel polynomials are not orthogonal 
in the conventional sense, i.e. over a real interval (see section 3.6 below for details). 

Some of the properties of the fourth polynomials have been studied in the current mathematical 
literature such as Refs. [2S],[2Z]- [30] ■ Their weight function is calculated from Eq. (|2.10p as 

w'-P^'^^ix) = (a;2 + i)-Pe9'^""'^ . (2.15) 

This weight function has first been reported by Routh |,24j . and independently Romanovski [25j . 
Notice that for a vanishing q, uj'^P'''\x) becomes equal to the student's t distribution [46] and acts 
as that very probability distribution. The a;^^'°-'(a:) case stands for the popular Cauchy (or, Breit- 
Wigner) distribution. Within the light of this discussion, the weight function of the Romanovski 
polynomials has been interpreted in Ref. [IB] as the most natural extension of the student's t 
distribution. 

The polynomials associated with Eq. (|2.15p are called after Romanovski and will be denoted by 
r'£'''\x). They have non-trivial orthogonality properties over the infinite interval x €z [—00, +00]. 
Indeed, as long as the weight function decreases as x~'^p, hence integrals of the type 

/+00 
w'^P'"^ (x)i?(P'«) {x)Rl';f^ {x)dx, (2.16) 
-00 

are convergent only if 

m + m' <2p-l, (2.17) 
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meaning that only a finite number of Romanovski polynomials are orthogonal. This is the reason 
for which one speaks of "finite" Romanovski polynomials. The orthogonality theorem has been 
proved in Refs. [571 ISZ] ■ The differential equation satisfied by the Romanovski polynomials reads 
as 

(1 + ^ ' + (2(-p + l)x + q) " ^ ' - {n{n - 1) + 2n(l - p))i?(f'«Hx) = . (2.18) 

a^x dx 

In the next section we shall show that the Schrodinger equation with the hyperbolic Scarf potential 
reduces precisely to the very Eq. (|2.18p . 

2.3 Bochner's classification. Explicit polynomial construc- 
tion. 

The earliest classification of the solutions to the generalized hypergeometric equation is due to 
Bochner [?7] (see also Refs. [341 1571 155] for more recent works) and based on the form of a{x) 
which can be rephrased in terms of the cr(x) roots. The presentation in this section closely follows 
Refs. [3ni|37]. Within this scheme one finds the following five polynomial classes : 

1. 'j{x) is a constant: 

The canonical form of the generalized hypergeometric equation is 

H"{x) - 2axH'{x) + \H{x) = 0, (2.19) 

where a is a parameter. The solutions are the generalized Hermite polynomials {iJ^(a;)} {a — \ 
characterizes the ordinary Hermite polynomials), and their weight function is 

J'^\x) ^ e-°"'\ (2.20) 
The case a — 1 provides the standard orthogonality relation. 



oo 



e-'-' Hn{x)H^,{x)dx = 5nn', V e {0,1,2,...} . (2.21) 



The lowest Hermite polynomials are then obtained as: 

Ho{x) - 1, 
Hi{x) = 2a;, 
H2{x) = -2 + 4a;2, 



Hz{x 
Hi{x 



= -12a; + 8a;^ 

„2 



= 12 - 48x^ + IGx^ (2.22) 



and are represented in Fig. 12.11 

2. <t{x) is of first degree: 

The canonical form of the corresponding hypergeometric equation is obtained as 

xL"{x) + yi{x)L'{x) + \L{x) = 0. (2.23) 

For 2/1 (x) = —ax + /3 + 1, a = 1, and /3 being a real number, the latter equation coincides with 
the equation for the associated Laguerre polynomials, |L^i^'''''(a:)|, while for a — \, and /? = it 
describes the ordinary Laguerre polynomials. The weight function is given by 

w("''^)(a;) =x^e-"^. (2.24) 
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f (X) 




Figure 2.1: Hermite polynomials. 

f (X) 




Figure 2.2: Laguerre polynomials for /3 = 1. 



For a, f3 > the orthogonality integral within the interval x S [0, oo) reads 

x^'e-''^L(„"^^\x)L^^'''\x)dx^6nn', V n,n'e {0,1,2,...}. (2.25) 

From now onwards we shall suppress the first upper index in L'^^'^^x) for simplicity. The lowest 
five associated Laguerre polynomials for, say, /3 = 1, are now calculated as 



L«(.x 



1, 

2-x, 

-{6 — 6x + x^), 

-(2A~36x + 12x^ - x^), 
6 

1 



— (120 - 240a; + 120a;2 - 20x^ + x^). 
24 



They are shown in Fig. 12.21 

3. cr^x) is of the second degree, with two different real roots: 

The canonical form of the corresponding hypergeometric equation is obtained as 

(1 - x^)P"{x) + yi{x)P'{x) + \P{x) = 0. 



(2.26) 



(2.27) 



The latter equation coincides with the Jacobi one for yi (a;) = a — /3 — (a + /3 + 2)a; where a, /3 G M 
are the polynomial parameters. The Jacobi polynomials, |pi"'''''(x)|, are then defined by the 
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f (X) 




4 



X 



-1 



1 



Figure 2.3: Jacobi polynomials for the parameters a = 1 and (3 — 2. 



Rodrigues formula in terms of the following weight function: 



(2.28) 



Furthermore, one has to restrict the parameters to a, /3 > —1 in order to ensure orthogonality in 
the interval x S [—1, +1] according to 



Several particular cases have received their proper names. These are: 

• Gegenbauer when a — [3^ 

• Chebyshev I and II when a — (3 — ±1/2, 

• Legendre when a = (3 — 0. 

The lowest five Jacobi polynomials for the toy values a = 1, and (3 = 2 ol the parameters are 
explicitly calculated as: 




(2.29) 




1, 

i(-l + 5x), 

3 + 9(-1 + x) + ^(-1 + .t)2, 

4 + 21(-1 + x) + 28(-l + xf + y (-1 + x)^ 

165 

5 + 40(-l + x) + 90(-l + xf + 75(-l + xf + —(-1 + xf. 



(2.30) 



They are shown in Fig. 12.31 



The lowest Gegenbauer polynomials and for a = 1 read: 



C^^\x)^l, 
C^f\x) = 2a;, 
Cf^(2:) = -1 + 4.t2, 
&^\x) = -4a; + 8a;^ 
&f\x) = -1 - 12a;2 + 16a;^ 



(2.31) 
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Figure 2.4: Gegenbauer polynomials for a — 

f (X) 
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Figure 2.5: Chebyshev I polynomials. 

They are graphically displayed in Fig. 12.41 
The lowest Chebyshev I polynomials read: 

To{x) = 1, 

Ti{x) — X, 

T2{x) ^ -I + 2x^ , 

T^ix) = -3x + Ax^, 

T4{x) = -1 - 8x^ + Sx'^ . 

They are shown in Fig. 12.51 
Next, the lowest Chebyshev II polynomials are: 

Uo{x) = 1, 
Ui{x) = 2x, 
U2{x) = -l + Ax^, 
U3{x) = -4a; + 8a;^, 
Ui{x) -1 - 12x2 ^Q^4^ 

They are graphically displayed in Fig. 12.61 
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Figure 2.7: Legendre polynomials. 



Finally, one finds the lowest Legendre polynomials as 



Po{x 

P2{X 

Pi{x 



1, 

X, 



2(1-3^')' 

-i(3-5x3), 



They are shown in Fig. 12.71 



(2.34) 



Related to the Legendre polynomials are the associated Legendre functions which are defined as 



p,r(2:) = (-i)™(i-^')^ 



dx' 



:Pn{x). 



(2.35) 



4. (j{x) is of the second degree, with one double real root: 

The canonical form of the corresponding differential equation is obtained as 

x^y"{x) + yi{x)y'{x) + Xy{x) = 0. 
For yi{x) = ax + /?, with a, /3 e M, the weight function is given by 



(2.36) 
(2.37) 
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Figure 2.8: Bessel polynomials for a — P — 2. 



The polynomials |yrr'^''(a;)| with a = (3 = 2, are called after Bessel [48]. They are not orthogonal 
in the conventional sense within a real interval. By considering x € C these polynomials have 
special characteristics. Only when the integration contour is the unit circle in the complex plane 
can one design orthogonality integrals and with respect to the "weight function" . To see this 
notice that Eq. (|2.36p with a = f3 = 2 can be rewritten to give 

{x^e-iyl,{x)y = n{n + l)e-* y„(a;) . (2.38) 
An integration by parts counter-clockwise around the unit circle leads to 

n{n + l) / ym{x)yn{x)e^idx = / {x'^e^iy'^{x))'yr,t{x)dx ^ - / x^e~^ y'^{x)y'^^{x)dx. (2.39) 
Ju Ju Ju 

Upon interchanging m against n and a subsequent subtraction, one arrives at the orthogonality 
relation 

(2.40) 

Ju 

The lowest Bessel polynomials are 





/ 2/m 


{x)yn{x)e ^dx = Q. 




JU 




are: 






2/0 (a;) 


= 1, 




yi{x) 


= 14 


-X, 


2/2(2;) 


= 14 


- 3a; + 3a;^, 


2/3(2;) 


= 14 


- 6a; + 15x^ + 15a;^, 


2/4(2;) 


= 14 


- lOx + 45x2 + 105a;^ + lOSx"* 



(2.41) 

They are displayed in Fig. 12.81 

5. cr(x) is of the second degree, with two complex roots: 

This is the case of prime interest to the thesis. The generalized hypergeometric equation for 
this choice of (j{x) takes the form 

(1 + x'^)R"{x) + yi{x)R'{x) + XR{x) = . (2.42) 

Taking yi{x) as yi{x) — (2/5 + l)a; + a, and a, /3 e M, we encounter a family of polynomials which 
we denote by |i?l'''"'*(a;)|. These are the Romanovski polynomials. The corresponding weight 
function (also mentioned before in Eq. (|2.15p ) reads 

Oj(l^^°'\x) = (1 + a;2)/5-ig-ata„-i(x)^ (2.43) 
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Figure 2.9: Romanovski polynomials for a — 1,(3 — — 1. 



Notion 


Symbol 


Weight Fimction: uo{x) 


Density: a{x) 


Interval 


Parameters 


Generalized Hermite 






1 


(— oo, oo) 




Generalized Laguerre 






X 


[0,oo) 




Jacobi 




(1 - x)°'{l + xf 


(l-x^) 


[-1,1] 


(«,/?> -1) 


Legendre 


Pn{x) 


1 


{1-x') 


[-1,1] 




Gegenbauer 


Cn (x) 


(l-a;2)A-i 


il-x') 


[-1,1] 


(A > -1/2) 


Chebyshev I 


Tnix) 




{1-x') 


[-1,1] 




Chebyshev II 


Unix) 




[l^x') 


[-1,1] 




Generalized Bessel 


V'n^'\x) 


X e ^ 


x^ 


complex 




Romanovski 


R't"\x) 


(1 + 2.2^/3-5g-atan~^(a;) 


(l + x^) 


(— oo, oo) 


-(] > n 



Table 2.1: Comparison among the different families of polynomials. 



For P — —p +1/2, and a — —q Eq. (|2.43p coincides with Eq. (j2.15p . The latter polynomials have 
special orthogonal properties and are studied in greater detail in the next chapter. 
The lowest Romanovski polynomials and for the toy values (3 — —1, and a ~ 1 are: 

4-'^'\x) = i, 

r["^'^\x) = -l-9a:, 
Ri'^'^\x) = -6 + 16a; + 56a;^ 
Ri'^"^\x) = 16 + 84a; - 126a;2 - 210a;^ 

R[~^'^^ (x) = 20 - 240a; - SeOa;^ + 480a;^ + 360a;^. (2.44) 

Some of them are shown in Fig. 12.91 In summary, the polynomial solutions of the generalized 
hypergeometric equation that permit for a Rodrigues representation fall into five classes. Three 
of them correspond to the celebrated classical orthogonal polynomials of Jacobi, Laguerre, and 
Hermite (recall that Gegenbauer, Chebyshev I and II, and Legendre appeared as special cases of 
the Jacobi polynomials). The principal characteristics of those five classes are collected in Table 
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Chapter 3 

Orthogonal polynomials in 
quantum mechanics. 



As already mentioned in the introduction, the classical orthogonal polynomials shape the exact 
solutions of a variety of quantum mechanical potentials. This is so because the Schrodinger 
equation with one of those potentials can be transformed into Eq. (j2.ip upon an appropriate 
change of the variables. This section is devoted to the one-dimensional Schrodinger equation 
(which equally well can be the radial part of an S wave 3D problem). In the following, we shall 
assign each polynomial in Table 12.11 to a potential that is exactly solvable in terms of that very 
polynomial. 



3.1 The Schrodinger equation with a radial potential. 

Two particles within a potential that depends only on the relative distance, r, 

T/(ri,r2) = F(r), r^ln-raj, (3.1) 

with ri, and r2 being the respective coordinates of first and second particle, are described by a 
time-independent Schrodinger equation of the form of a single particle within the center of mass 
frame, 

-—W'^-f{r) + V{r)^{r) = E^{r), (3.2) 

where fi stands for the reduced mass. This equation can be solved by standard techniques like 
separation of the variables, in which case the solution factorizes into a radial function, u(r), 
and an angular part, Q(6)^{ip), which for central potentials is given by the spherical harmonics, 
&{eMip)^YhniO,ip), i.e. 

*(r) =u(r)y,™(0,(^). (3.3) 

In substituting the latter wave function in the three-dimensional Schrodinger equation, amounts 
to the well known one-dimensional radial equation 

fi^ 1 d f 2 d . .\ 1(1 + 1) , . ^ , . 

-2il^d-r [' Tr''^'^ ) + 2^1 + ^^^^^^'^^ =^ ^^'^^ 

In making the coordinate and the wave functions dimensionless upon rescaling, 

^=d' "W = ^' (3.5) 
where d is an appropriate length scale, Eq. p.4p takes the form 
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Here, 

v{z) = V{r)l{h^/2n(P), e = E/{h'^/2fi(f), (3.7) 
and the normalization is defined by the integral 



oo 

2 



0>(z)V'*(z)rfz = 1. (3.8) 



In terms of the new variable, R{z), introduced via 



i^{z) ^ (3.9) 



Eq. (|3.6p becomes 



33ji?(z)+ (^--^+v(z)-e)i?(z) = 0. (3.10) 



With that the normalization integral changes to 

/>oo 

/ R{z)R*{z)dz = \. 
Jo 



(3.11) 



Casting the three-dimensional Schrodinger equation into the form of Eq. p.lOp . brings two im- 
portant advantages. The first one is that all variables are dimensionless, which allows for an 
easier mathematical treatment. The second advantage is that Eq. p.lO[) takes the form of ID 
Schrodinger equation with an effective potential defined as 

vu{z)=viz)+^-^^^. (3.12) 

Finally, not to forget the -R(O) = boundary condition, telling that the wave function has to vanish 
at the origin if it is not to penetrate the centrifugal barrier. In this fashion, the 3D problem has 
been replaced by an equivalent ID problem. 



3.2 Orthogonal polynomials in exact wave functions. 

The Schrodinger wave functions, ijj{z), in diagonalizing an Hermitian differential operator, are 
known to constitute a complete orthogonal set. One can try to reduce the radial Schrodinger 
equation (written in the variable z) to an appropriate polynomial equation (written in the variable 
x) by means of the following change of variables: 

z = f{x), R{f{x))^g{x), g{x) ^ JJ^^) F^^^-^\x) (3.13) 

V dx 

If this substitution turns out to be successful, the differential equation for Fi°''^\x) will be a 
version of the generalized hypergeometric equation whose solutions have been presented in the 
previous chapter. When expressed in terms of polynomials, the orthogonality of the Schrodinger 
wave functions translates into orthogonality between the involved polynomials according to 

r+oa A2 

/ Rn{z)Rn'{z)dz = / g„{x)gn'{x)df{x) 
JO Jh 

= J''^l^\x)F^^^'^\x)Fl^'^\x)dx. (3.14) 

The parameters a, and /3 are some functions of the potential parameters. Occasionally, the 
price of casting the Schrodinger equation in the form of Eq. (|2.5p is that a and (3 acquire an 
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n-dependence. In such cases, as we shall see below, Eq. p.l3p may not hold valid and the orthog- 
onality of the wave functions may not amount to the orthogonality integral of the polynomials 
with free parameters. In the following we employ the idea of Eq. (|3.13p in order to solve the 
Schrodinger equations for a variety of potentials. We show that the Hermite polynomials shape 
the solutions of the ID oscillator, the Laguerre polynomials define the wave functions of a par- 
ticle within the 3-D-oscillator- and the Coulomb wells. The Jacobi polynomials solve exactly the 
hyperbolic Rosen-Morse (Rosen-Morse II) potential, while the e"'' barrier is solved in terms of 
the Bessel polynomials. The results listed in Table 1.1 above have been obtained pursuing same 
path. 



3.3 Hermite polynomials and the harmonic oscillator. 

The simple harmonic oscillator in quantum mechanics is a very well known example elaborated in 
all the standard textbooks. Its solutions are given in terms of the Hermite polynomials. In order 
to see this we write down the ID Schrodinger equation, which reads: 

-^-^^iy) + Viymy)^Eq^{y). (3.15) 

Here, m is the particle's mass, V{x) stands for the potential, and E is the energy of the system. 
The simple harmonic oscillation potential is given by 

Vnoiy) = IrncuV: (3.16) 

where lo is the frequency of the oscillations in classical mechanics. Substitution of the latter 
equation into Eq. p.lSp gives 

^cj^iy) + ^(e- lmc.v) Hy) = 0. (3.17) 



dy^^^"" h? \ 2 
The two parameters of the problem (to, lu) provide a unit of length for the problem: 



^, . = Hy = zd) = ^. (3.18) 



In performing a substitution of the form in Eq. (|3.13D and given by 



z^f{x)^x, ^{z^x)^g{x), g{x)^NnV^H{x) , (3.19) 

where iV„ is a normalization constant, a new equation is obtained (after substitution and reorder- 
ing): ^ 

±^g{x) + {e-x^)g{x)^0, (3.20) 
where e = E now is dimensionless one finds 

H"{x) - 2xH'{x) + (e - l)H{x) = 0, (3.21) 

which is the Hermite equation 

H"{x) - 2xH'{x) + \H{x) = 0, (3.22) 

with A G M a constant to be determined. This equation admits a polynomial solution of degree n, 
Hn{x), only if A — 2n. Since we are interested in such solutions, we conclude 

e„-l = 2n, ne {0,1,2,...}, (3.23) 
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Figure 3.1: The simple ID harmonic osciUator, V{x) 
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^x^. Here muj = 1. 



I.e. 



^2 1 1 



(3.24) 



The orthogonahty integral between the wave functions amounts to the orthogonality integral of 
the Hermite polynomials 



ipn{z)ijM'iz)dz 



(3.25) 



The harmonic oscillator potential finds application in the description of vibrational modes in 
nuclei, atoms, molecules, and crystal lattices. This potential is shown in Figure [ST] 



3.4 Laguerre polynomials in 3D oscillator and Coulomb wells. 
3.4.1 The 3D oscillator. 

Harmonic oscillations in three dimensions lead to a differential equation whose solutions are the 
associated Laguerre polynomials. One way to solve this problem is to consider the three dimen- 
sional oscillation as the result of oscillators in x, y, and z directions. The solutions can then be 
expressed in terms of the product of three Hermite polynomials. However, given the rotational 
symmetry of the problem, the most natural coordinate choice are the spherical coordinates in 
which case the complete solution factorizes in a radial function and spherical harmonics in accord 




Figure 3.2: Three-dimensional harmonic oscillator potential (dashed line), V{r) 



2J1 



Effec- 



tive potential (solid line), Veff{r) 
moj = 1, and 1 = 1. 



2^2 _|_ 1(1+1) includes the centrifugal barrier. Here 
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with Eq. (|3.3p . The 3D oscillator potential reads 

V{r) = ^rriLU^r'^, (3.26) 

and is shown in Fig. 13.21 Changing variables to z — r/d and substituting into Eq. p.6p results in 

(f 2 d 1(1 + I) A , , , 

-- +^^+z^\^p{z)^eyj{z). (3.27) 
dz^ z dz z'^ J 

The latter equation has to have well behaved asymptotic solutions at the origin and infinity. For 
z — » oo it reduces to 

d2 



-— +z^jV^(z) = eV'(z), (3.28) 

meaning that ij] k, would be a good approximation. Next, we inspect the behavior of 

Eq. (|3.27p at the origin, z ^ 0, where it reduces to 



d^ 2d + 
dz"^ z dz z^ 



i){z) = 0. (3.29) 



Taking '0(-2^) = z* as a test function, we find that it would be a solution for t = I, or, t = —I — 1. 
The former value is the correct one since it is finite at the origin. Now we shall find the exact 
expression for the solution applying the standard techniques we used before. 
The above asymptotic behaviors suggest to try as an ansatz 

z^{z) = R{z) = K„z^+hr^F{z), (3.30) 
where Kn is a normalization constant. Insertion of the latter expression into Eq. p.27p leads to 

(-1 + 2z2 - 2;9)i^'(z) - zF"{z) + i ^— + z{2 + 2/3 - e) j F(z) = 0, (3.31) 

where /? is a free parameter of a choice suitable for simplifying the resulting equation. Next we 
notice that in terms of the new variable, x, introduced as z = fix) — \/x, the ansatz in Eq. (|3.30p 
can be converted to R{^/x) = g{x) = Nnx'i^ie'^ L{x) with F{z = y^) = L{x), and Nn a new 
normalization constant. This expression is of the type in Eq. (j3.13p in so far as 



g{x) = NnVxf^e-^L{x)/^ -^V^, 



(3.32) 



where one recognizes the square root of the weight function of the associated Laguerre polynomials 
in front of L(x). The difference between iV„ and Kn accounts for possible constants emerging from 
the inverse of the derivative of f{x). In effect, one arrives at 

xL"{x) + (/3 + 1 - x)L'{x) + [^^l^—^^l±^l^ + £ _ ^ _ 1^ L{x) ^ 0. (3.33) 

Now we make use of the freedom in (3 to nullify the singularity by setting (3 = 1 + ^ and thus, we 
are left with the equation for the associated Laguerre polynomials, equation (|2.23p . Therefore 

4/32 - 4/(Z + 1) - 1 = 0, (3.34) 

allows to identify Eq. p.33p with the equation for the associated Laguerre polynomials, 

xLl{x) + (Z + i + 1 - x)L'n{x) + nLn{x) = 0, (3.35) 



24 



whose solutions are Ln{x) — L\ ' (x). In effect, the radial part of the Schrodinger wave function 
for the 3D oscillator is given as 



Back to the z variable, the Schrodinger wave function is now obtained in its final form as 



(3.36) 



(3.37) 



The orthogonality integral between the wave functions recovers the orthogonality between the 
Laguerre polynomials according to 

Rn{z)Rn{z)dz = r KnKn>Vx'+h--Li!'^^\x)Vx'+-2e'-4':^~'\x)dx = S^n' , (3.38) 



'0 "'0 

where dz = d^/x. The energies are now given by e„ = 2(2n + / + §). The reason for which 
the associated, and not the ordinary Laguerre polynomials appeared in the solution of the 3D 
oscillator is that the angular momentum, I, requires a polynomial parameter. 



3.4.2 The hydrogen atom. 

In this subsection we study the hydrogen atom, which is a two-body system consisting of one 
proton and one electron. The interaction is governed by the Coulomb potential. 



V{r) = - — , Z = l, 
r 



(3.39) 

as displayed in Fig. 13.31 Substituting the latter equation for the potential in Eq. (j3.10p leads to 

- ^R{z) + {^^^^^ ^^^^ = °' ^ " 

where the length scale d, has been chosen as the Bohr radius, d = /i^/(/ie^), and e stands for 
e = 2Ed/e^. We are going to solve this equation by means of an appropriate variable substitution. 
In order to find it we first study the asymptotic behavior of the solutions at both origin and 
infinity. For finite z the latter equation is equivalently rewritten as 



^^^^^^ + (^^^T^ - 2 - ez ) Riz) = 0, 



(3.41) 



V(r) 




Figure 3.3: Coulomb potential (dashed line), V{r) — —^f-- The effective potential (solid line), 
Veff{r) = — h '^'^"'"•^ includes the centrifugal barrier. Here Ze^ — 10 and I — 1. 
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while near the origin, where i?(0) = 0, and 1/z << 1/z^ one finds 

1(1 + 1) 



dz^ 



R{z) w 0. (3.42) 



The latter equation can be solved by a function of the type, R{z) = z*. This would restrict t to 
t = I + 1, and t = —I, respectively. The second solution is not acceptable because it is not finite 
at the origin, and we are left with the first one. Next we study the z —> oo asymptotic solution. 
In this case Eq. p.40p becomes 

As long as we here are interested in the bound states, e < 0, which amounts to 

d^ 

^i?(z) - \e\Riz) = 0. (3.44) 

The latter equation is solved by R{z) — CeV^^ and R{z) — Ce~V^^. As long as we want the 
wave function to go to zero for large z, the first solution has to be dismissed. Next one changes 
variables in Eq. p.40p according to 



z^f{x) = -, K^2^\, i?(-)=.g(x), (3.45) 



K \K 

where the factor 2 has been taken for convenience. This yields 



d^ 1(1 + 1) 1 1\ , , , , 

+ + - - I 9{x) = 0. (3.46) 



Taking into account the correct asymptotic behavior revealed above, the solution of the latter 
equation can be assumed as g{x) — x'-^^e^^ L{x). In result, one arrives at 

xL"{x) + {21 + 1 + 1- x)L'{x) + L{x) [ ^= - / - 1 I = , (3.47) 



which coincides with the equation (|3.35p for the associated Laguerre polynomials provided a — 1, 
/3 = 2Z + 1, and \ = n. With this in mind it is easy to verify that 

[n + l + 1)'' 

Recalling the definition of (e = 2Ed/e'^), we obtain the spectrum as 

- - 2nHn + i + i)^ - ^'-''^ 

At that stage one defines iV = n + Z + 1 as the principal quantum number, and realizes that while n 
can be any non-negative integer, the angular momentum is restricted to a finite number of values 
according to ^ = — 1 — n, 

/ = 0,1,2,...,A-1. (3.50) 
In effect, one finds the well known degeneracy patterns in the spectrum of the hydrogen atom, 

4 

Eni — ~ 2K^n'^ ' meaning that k depends on n. The solution for the radial equation (|3.47p is then 
concluded as 

^ = V„;(^) = C,ux'e-^-L^^'+^\{x), (3.51) 

X 

where use has been made of Eq. (j3.9p . The latter expression equivalently rewrites to 



i^ni{x) = C„iVa;2'+ie--42'+i)(x)/Vi. (3.52) 
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Therefore, -0„;(a;), is given in terms of the associated Laguerre polynomials, with 
Cni being the normalization constant. However, the wave function under consideration is not 
of the form in Eq. (j3.13p because y/x ^ (x/k)' . The consequence will be that the orthogonality 
integral of the wave functions 

I / w 

aiC„,,y^Kz)2'+ie--"-42'+i)(^„z)y^K,z)2'+ie--.-4^'+i)(«:„,z)z2^ = <5„„, , (3.53) 

does not coincide with the orthogonality integral of the Laguerre polynomials with 
free parameters because it contains the additional first power of z. It is important to be aware 
of the fact Eq. (|3.53p describes orthogonality between states bound within different potentials 
corresponding to different factors, ^-^^ — j: versus ^-^= — j:, of 1/x in Eq. (j3.46p . In the textbook 

on Mathematical Methods in Physics by Arfken and Weber (second reference in [31]) this phe- 
nomenon has been also attributed to the dependence of k on the degree of the polynomial via the 
energy (see Exercise 13.2.11 there). As we shall see below, such a behavior is much more general 
and can occur also when x is neat but the parameters carry an n dependence. 

3.5 Jacobi polynomials in Rosen-Morse II. 

The Jacobi polynomials appear in several physics problems ranging from classical electrodynamics 
to quantum mechanics. We here focus on the hyperbolic Rosen-Morse potential given by 

52 

v(z) = a^ + — - a(a+l) sech^(z) + 26 tanh(z), (3.54) 



and displayed in Fig. 13.41 The corresponding ID Schrodinger equation is 

(PRiz) / 52 \ 

, \ ' + + aia + 1)(1 - tanh2(z)) - 26tanh(z) + e Riz) ^ 0, (3.55) 

dz-^ \ / 

with e = a? — t, and s.ec\\^(z) = 1 — tanh2(z). Changing variable to x = tanh(z) meaning z — f{x) 
with f{x) = tanh^^(x), results in 

, ^.dPqix) „ dq{x) f , , 2bx + e\ , , „ ,„ 




Figure 3.4: The Rosen Morse II (hyperbolic) potential, V{r) = + — a(a + l)sech2(r) -|- 
2&tanh(r). Here a = 10 and 6 = 10. 
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with g{x) being defined as g{x) — i?(tanh ^x). By means of the substitution, g(x) ~ (1 + a;) 2 (1 
x)^P[x), the latter equation becomes 

(l-x^) , ; ^ +(/3-a-a:(2 + a + /?)) ^ ' ' 



dx^ dx 

(3.57) 



where a, (3 are free parameters to be used to simplify the equation. Specifically, one makes use of 
the freedom in a, (3 to nullify the singular term which restricts the parameters to: 

2fo-Y + Y=0, (3.58) 

^---f-e = 0. (3.5.,) 

Next one requires the constant multiplying P{x) to be of the form 

a(a + 1) - ^ii±^ _^_^_^^A„ = n(l + n + a + /3). (3.60) 



The latter equations are resolved b 



P — a~n~\ = fJ-n, (3.61) 



n — a 



a = a — n = z/„, (3.62) 

n — a 

and 

e„ = ^-(a-n) (3.63) 

[a — ny 

With that, equation p.57p can be identified with the Jacobi form of the generalized hypergeometric 
equation. 



(1 - x^)^-^ + (a - /3 - x(2 + a + /3)) — +n{n + a + (3+ l)P{x) = 0, (3.64) 
da;^ ax 



whose solutions are the Jacobi polynomials [pL''''\x)}. In effect, the hyperbolic Rosen-Morse 
potential is solved exactly by 

g„(x) = ^{l~xy^{l + xY^Pi'^-'^-\x), (3.65) 

and in accord with Table 1.1 (when translated to our notations). This wave function is not of the 
form in Eq. (j3.13p . As a consequence, the orthogonality integral between the wave functions does 
not recover the orthogonality between the Jacobi polynomials with free parameters 

as visible from 

+00 



7V„iV„.i?„(z)i?„.(z)dz - j V(l-2:)^"(l + x)-"P^^"'''")(x) 

v/(l - a;)^"'(l + x)^"'P^^''''"'^(x)dtanh"^(x) . (3.66) 

The culprits for this are the n dependent polynomial parameters. This is not to remain the only 
example for such an anomalous behavior. 
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3.6 The Bessel polynomials in spherical waves phenomena. 



The Bessel polynomials are special in the sense that their orthogonality is achieved by integration 
over a contour in the complex plane. In the context of quantum mechanics an infinite barrier 
(Fig. 13. 5p of the type Ae""^ {A, a being constants) will lead to Bessel's differential equation. Using 
this barrier in Eq. (|3.10p for the case I = results into 



Changing variables to a; = e~"^ leads to 



+ (e + Ae"")i?(z) = 0. 



A substitution of the type g{x) — e '^x^y^x) leads to 



~^ dx 




e 













g{x) = 0, 



^(fy{x)_ 
dx^ 



+ 2{x + l) 



dyjx) 
dx 



1 - 4e 4-46 



4a;2 



yix) = 0. 



If the singular term is to vanish, then 



4 - 46 = 
4e 



n{n + 1). 



(3.67) 

(3.68) 
(3.69) 

(3.70) 

(3.71) 
(3.72) 



From the above conditions we find 6 
differential equation: 

x2 



,d^y{x) 



1, e — J + n{n + 1) and Eq. (|3.70p becomes the Bessel 
dy{x) 



+ 2(a;+ 1)- 



n(n + l)y{x). 



(3.73) 



dx"^ ' dx 

Here, y{x) stand for the Bessel polynomials. Take notice that contrary to the previous examples the 
barrier potential is exactly solvable only when the parameters have been fixed by the condition 
^ = 1. This is an example for a potential that is exactly solvable when the parameters obey 
constraints. The Bessel polynomials have been brought to attention by Krall and Frink [15] and 
are related to the following wave equation (in spherical coordinates) 



1 / 52 

sin2(6l) 



1 



d 



1 d^u{r,9,ip,t) 



du{r, 9, ip,t) 
39 



(3.74) 
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Figure 3.5: A barrier V{r) ~ Ae'"' that leads to a Schrodinger equation with solutions in terms of 
Bessel polynomial's and for the toy values a — —10, and A=\. 



29 



If the latter equation is solved by separation of variables, the radial part /(r) is found to satisfy 
the following differential equation: 

r^fip. + + fcVVW = n{n + l)/(r). (3.75) 

dr-' dr 

Here kc = oj, and e*"^* is a plane wave. For r — ks, Eq. p.75p becomes 



s'^^ + 2s^ + s'm^nin + l)m, /(s) = / (r = . (3.76) 



The latter equation can be transformed into Bessel's equation by the substitution /(s) = 52 J(s), 
which yields: 

^2 cPj(5) ^ ^dJjsl _^ ^2 j(^) ^ ^ 1/2)V(.). (3.77) 
as^ ds 

For integer n, the solutions of Eq. (|3.77p are the Bessel functions of half-integral order and are well 
known. However, Eq. p.76p can also be subjected to the transformation /(s) = w{s)/s, giving 

s' (^^J^ + ^(*)) - + (3-78) 

Upon introducing the variable, z = ikr = is, imaginary and real parts of the solution represent 
traveling waves. Because of that we admit in Eq. p.78p z = is = ikr, and w(s) = e~^y{z) = 
e~^^y{is), and obtain 

.^(^-2^).n(n+l),(.). (3.79) 

For integer values of n, Eq. p.79p has solutions which are polynomials in 1/z. Therefore the final 
substitution should be, x — 1/z = 1/ikr, which allows to obtain the differential equation for the 
Bessel polynomials {Unix)}, 

x^^^ + (2x + 2)^^=n(n + l)y„(x). (3.80) 
dx"^ X 

The full solution of the spherical wave equation is then given by 

u{r, e, If, t) = r-^P^'{cos{9)) sm{frvfi - a)e'^'^'-'''''>yn{l/ikr), (3.81) 

where yn{x) = yni^/ikr) is a Bessel polynomial, and kc = lo. The real and imaginary parts 
of p.8ip describe waves traveling in the radial direction with velocity c. In conclusion, spherical 
waves are equivalently described either in terms of a class of polynomials orthogonal over the unit 
circle, or in the standard way in terms of Bessel functions. 



3.7 Romanovski polynomials in Rosen-Morse I. 

As already mentioned in the introduction, the Table 1.1 contains one more potential whose exact 
solutions require the Romanovski polynomials and this is the trigonometric Rosen- Morse potential 
(Rosen- Morse I). This case has been considered in great detail in Refs. [l9],[50j and will not be 
repeated here. Instead, I would prefer to briefly review the main properties of its solutions within 
the context of its relevance in quark physics (the presentation on this section closely follows 
Ref. [37]). This aspect appears especially important to me because the Romanovski polynomials 
were found for the first time in the solutions of that very potential while searching to construct a 
quark model that matches reality on nucleon excitations as part of the research project "Dynamics 
of baryon resonances" run by our group. The subject of the present thesis continues the study of 
the Romanovski polynomials started in Refs. [121 [SU] and extends knowledge on their properties 
by various new observations. 

The great appeal of Rosen-Morse I is that 
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besides the Coulomb potential, the trigonometric Rosen-Morse potential is to the best of 
our knowledge the only exactly solvable potential that relates to a fundamental massless 
gauge theory. 

Recall that the Coulomb potential is no more but the image in coordinate 
space of the propagator, — of the photon, the U(l) gauge boson of elec- 
trodynamics, as it appears in elastic scattering of charged particles. In 
a similar way, as it will be argued below, Rosen-Morse 1 can be viewed as 
the image in coordinate space of the propagator of the gluons , the SU{3)c 
gauge fields of the fundamental field theory of strong interactions, the 
Quantum ChromoDynamics (QCD) , as it appears in elastic scattering of quarks. 
Within this context, the Romanovski polynomials acquire the special status 
of major ingredients of the wave functions of bound quarks. 

The quarks, as is well known, are the constituents of strongly interacting particles, baryons, and 
mesons. In what follows, the presentation will be focused on the baryons constituted by the so 
called light flavors, u, d, and s, as are the nucleon {N), the A, the A, and their resonant excitations 
[51j . The nucleon is understood as a particle which can exist in two different states distinguished 
by their electric charges, the proton, p{uud), and the neutron, n(udd), and is said to be a charge- 
doublet. The A stands for a particle that can exist in four different states distinguished by their 
charges which are A'*""'"(uuu), A^{uud), A'^{udd), and A~{ddd) and is termed to as a charge- 
quadruplet. The A{uds) particle is neutral and a charge singlet. There are also charge triplets 
like S represented by S+(wMs),E'^(u(is), and E~(d(is), and one more charge-doublet given by the 
{S^ (dss), S'^(mss)} family. Compared to the N-, A-, and A- charge-multiplets, the E and S are 
less known and will be left out of consideration in the following. 

The problem which one is facing with the baryon resonances is the lack of an adequate sys- 
tematics and the resulting deficits as the prediction of a large number of unobserved states. In 
the standard quark models, different charge-multiplets are supposed to join to bigger families, the 
so called SU{6)sf x 0{3)l super-multiplets as displayed in Fig. 13.61 This figure reveals a strong 
overlap between the different super-multiplets and an apparent lack of degeneracy between the 
states belonging to same multiplet thus questioning the adequacy of the underlying classification 
scheme. In Refs. [Sl| the SU{6)sf x 0(3)l baryon classification scheme has been given up, the 
super-multiplets have been decomposed into N-, A-, and A-spectra and have then have been 
studied separately. 

The result was that to a surprisingly good accuracy, the nucleon excitation levels carry the same 
degeneracies as the levels of the electron with spin in the hydrogen atom, though the splittings of 
the former are quite different from those of the latter. Namely, compared to the hydrogen atom, 
the baryon level splittings contain in addition to the Balmer term also its inverse but of opposite 
sign. Same patterns are repeated by the excitation spectrum of the A(1232) particle, the most 
important baryon after the nucleon (see Figs. I3.7l3.8p . In this way baryons have been classified 
according to SU{2)j x 0(4) with / standing for isospin (a number that encodes the dimensionality, 
D, of the charge multiplets as £> = 2/+ 1). The appeal of the new classification scheme lies in the 
fact that no states drop out of the systematics, on the one side, and that the number of unobserved 
("missing" ) states predicted by it is significantly smaller than those of all preceding schemes. 
The observed degeneracies in the spectra of the light quark baryons have been attributed in 
Ref. [52] to the dominance of a quark-antiquark configuration in baryon structure. Within the 
light of these findings, the form of the potential in configuration space acquires crucial importance. 

In Refs. j4^>]S2il case was made that precisely the trigonometric Rosen-Morse po- 
tential provides the degeneracies and level splittings that are required by the light quark 
baryon spectra. 

The success of the trigonometric Rosen-Morse potential in quark physics is not accidental. It 
is due to the property of the latter to interpolate between the Coulomb- and the infinite well 
potentials. In order to understand this virtue one has to recall the QCD basics. 
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Indeed, the strong interaction between quarks within QCD is governed by interchangings of 
massless gauge bosons between quarks of three different "colors" , ("red", "blue", "green"). The 
above gauge theory predicts that the quark interactions proceed over one- or many gluon exchanges 
including gluon self- interactions (so called "non- Abelian" gauge theory) . The latter are believed to 
be responsible for the so-called quark confinement, where highly energetic quarks remain trapped 
but behave as (asymptotically) free at small distances. The QCD equations are nonlinear and 
complicated due to the gluonic self-interaction processes. Their solution requires employment 
of highly sophisticated techniques, such as discretization of space time, so-called lattice QCD. 
Lattice QCD calculations of the properties of hadrons associates the one-gluon exchange with 
a Coulomb like, ~ 1/r, potential, and predicts a linear confinement potential with increasing 
energy brought about by gluon-self interactions. The trigonometric Rosen- Morse potential in the 
parametrization of Ref. 53J where the a parameter has been identified with /, the relative quark- 
di-quark angular momentum (in units of fi,^ = 1 = 2/i, and d a suited length scale making the z 
variable dimensionless), 

Wti?M(^) = -2&cot(z) + /(/ + l)— 4—, (3.82) 

sm (z) a 

has precisely the properties required by lattice QCD. It captures the essential traits of the QCD 
quark-gluon dynamics in interpolating between the Coulomb potential (associated with the one- 
gluon exchange) and the infinite well potential (associated with the trapped but asymptotically 
free quarks) while passing through a linear confinement region (as predicted by lattice QCD) (see 
Fig. 13. 9p . In order to see this one has only to perform the Taylor expansion of the potential of 
interest, 

26 2& 1(1 + 1) 1(1 + 1) 2 

z o z J_o 

This expansion clearly reveals the proximity of the cot(z) term to the Coulomb- plus linear 
confinement potential, and the proximity of the csc^(z) term to the standard centrifugal barrier. 

In this sense, Rosen-Morse I can be viewed as the image of space-like gluon propagation 
in coordinate space. 

The great advantage of the trigonometric Rosen-Morse potential over the linear- plus Coulomb 
potentials is that while the latter is neither especially symmetric, nor exactly soluble, the former 
is both, it has the dynamical 0(4) symmetry (as the hydrogen atom) and is exactly soluble. The 
exact solutions of the, now three dimensional, Schrodinger equation with VtRAiiz) from Eq. (|3.82p 
have been constructed in [53] on the basis of the one-dimensional solutions found in [13] and read: 

Mcot-'x) = (l + x2)-^e-™^°*"(^)ci^*"'"'^'™)(a;), (3.84) 

with X = cot(z). The C polynomials from [491 are Romanovski polynomials but with running 
parameters attached to the degree of the polynomial. The notations of Ref. [49j translate to the 
present ones as: 

^(-(«+0,7Jtt)(^)^^(p.,«.)(^)^ <Z„ = pn = n + l + l, 71 = 0,1,2,... (3.85) 

n + I + 1 

The wave function in Eq. (j3.84p is not of the type in Eq. (j3.13p . a reason for which the orthogonality 
integral between the wave functions will not recover the orthogonality of the Romanovski 
polynomials vifith free parameters as visible from 



Mz)i'n'{z)dz = r 7l.(P"^?")(x)7V„i?i^'-«")(x)77«(P"''?"')(x)A^„-i?i^''«"'^(x)-4^ = 6n n' , 

(3.86) 

where the factor in the integrand comes from = — 1/(1 -f a;^). In this way, the Ro- 

manovski polynomials that enter the solutions of the Schrodinger equation with the trigonometric 
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Rosen-Morse potential seem to disobey the finite orthogonality prescription. A discussion of this 
behavior will be given below. Finally, the associated energy spectrum is found as 

e„ = (n + ; + l)2--— . (3.87) 

{n + l + 1)^ 

Therefore, the Romanovski polynomials have been shown in Refs. [ 491 153j to be important ingre- 
dients of the wave functions of quarks designed in accord with QCD quark-gluon dynamics. 
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7+ 5+ 3+ 7+ 1 3- 5- 7- 

Figure 3.6: Baryon resonances in the traditional quark model. Circles, bricks, and triangles 
stand for nucleon, A, and A states, respectively. Different colors mark different SU{Q)sf x 0(3)l 
multiplets. Notice the strong multiplet intertwining and the large mass separation inside the 
multiplets. (Courtesy M. Kirchbach) 
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Figure 3.7: The nucleon excitation spectrum below 2 GeV. (Courtesy M. Kirchbach). 
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Figure 3.8: The A excitation spectrum below 2 GeV. (Courtesy M. Kirchbach) 




Figure 3.9: The trigonometric Rosen-Morse potential (solid line) and its proximity to the 
Coulomb- plus lineal potential as predicted by lattice QCD (thin dashed line) for the toy val- 
ues ^ = 1, 6 = 50 of the parameters. 
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Chapter 4 



Romanovski polynomials in Scarf 
II. 



This chapter is devoted to the presentation of the original results obtained in the thesis. In the 
first section 4.1 I present the hyperbolic Scarf potential. The corresponding Schrodinger equation 
is solved in terms of the Romanovski polynomials and presented in subsection 4.1.1. The finite 
orthogonality of the Romanovski polynomials is discussed in section 4.2. Section 4.3 is devoted to 
the problem of an electron within a non-central potential, which is solved in spherical coordinates 
by separation of the variables. There, in sub-section 4.3.2 it is shown that the Romanovski 
polynomials solve exactly the polar angle equation and define new non-spherical angular functions. 
In subsection 4.3.3 a non-linear relationship between Romanovski polynomials and associated 
Legendre functions is established. Finally, in section 4.4 I solve the Klein-Gordon equation with 
scalar and vector potentials of same magnitudes and given by Scarf II. 

4.1 The hyperbolic Scarf potential. 

The hyperbolic Scarf potential can be viewed as the extension of the sech^(r) potential, or, bet- 
ter, of the original Poschl- Teller potential [Ij. Take notice that in the SUSYQM nomenclature 
(presented in Table 1.1) the names of Poschl and Teller are rather associated with the extended 
csch^(r) (here marked as Poschl- Teller 2). It seems that for the first time Scarf II has been con- 
structed in Ref. [5S]. It also has been encountered independently within the framework of the 
supersymmetric quantum mechanics [131 1221 123j while exploring the superpotential 

Umix) = a tanh(Q:a;) + b sech(ax). (4.1) 




Figure 4.1: The trigonometric Scarf potential (Scarf I) for the toy values of the parameters, a = 10, 
6 = 5, and a — 1. The horizontal lines represent the discrete levels. 
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Figure 4.2: The hyperbolic Scarf potential Vh{z) = + (6^ — — aa)sech'^{az) + b{2a + 
a)sech(az)tanh(az) with a = 10, 6 = 5 and a — I. Energy levels, e„ = ~ (a — na)^, are 
included. 

On the other side, it can equally well be approached from the perspective of the trigonometric 
Scarf potential, here denoted by Vt{z) and given by 

Vt{z) = -a^ + (a^ + fe^ - aa) sec^(az) - b{2a + a) tan(az) sec(az). (4.2) 

The exact solution of the Schrodinger equation with the trigonometric Scarf potential (displayed 
in Fig. l4.ip can be obtained along the line of the concepts of the previous chapter. It is well known 
and given in terms of the Jacobi polynomials, Pn^'^\z), as |13j, 

M^) = (l-sin(az))^(l + sin(az))^pi'°"°"^^''"°"""^^\z). (4.3) 

The corresponding energy spectrum is obtained as 

e„ = (a + an)^ - . (4.4) 

The trigonometric Scarf potential can be transformed into its hyperbolic partner, the so called 
hyperbolic Scarf potential, here denoted by Vh{z) and given by 

Vh{z) — + (b^ — — aQ;)sech'^(a2:) + b{2a + a)sech(az)tanh(Q!z). (4-5) 

Figure 14.21 visualizes the hyperbolic Scarf potential and its discrete spectrum. 

The Vh{z) potential has been obtained from Vt{z) in performing the following substitutions in 
Eq. g2D : 

a — > ia , a — *■ —ia , 

b^b. (4.6) 

a a 

Upon the above substitutions the energy changes to 

e„ = — (a — na)'^ . (4-7) 

In the following we shall show that n = 0, 1, 2, ... < a meaning that the number of bound states 
is finite. Yet, the most profound changes are suffered by the wave functions. Substitution of 
Eqs. (HH) into Eq. (g^]) results in 



ipni—i sinh(z)) = {1 + i sinh(z)) ^ (1 — i sinh(z)) 



I — i sinh(z) 
I + i sinh(2:) 



c„Pi^'+'^"^^^^"^''+""^V*sinh(^)) , (4.8) 



37 



where c„ is some state dependent complex phase, and where we took a = 1 for simpHcity. The 
latter expression can be cast into the form frequently mentioned in the literature [l3],[22], [36] . |45j . 

X — sinh(z) , p — a, a = b . (4-9) 

The latter equation gives the impression that the exact solutions of the hyperbolic Scarf potential 
rely upon Jacobi polynomials with complex indices and arguments. 

In this thesis the case is made that this needs not be so and that the above wave 
functions can be expressed in terms of the real Romanovski polynomials. 

4.1.1 The polynomial equation. 

The Schrodinger equation for the potential of interest when rewritten in a new variable, x, in- 
troduced via an appropriate point canonical transformation [56j . [57] . taken by us as z = f{x) = 
sinh~^ x, is obtained as: 

o.d'^gix) dg{x) f-9 + a{a + l) b{2a + l) \ ^ 

with g{x) = 'i/;(sinh^^(x)). Inspired by Eq. (|4.9p we now test the following substitution in Eq. (|4.10p 
g{x) = (l + a;2)4e-t'^""'(^)i:)('5^")(a;), x = sinh(z), (4.11) 
In effect, Eq. (j4.1ip reduces to the following equation for D'-^'°'\x), 

(1 + x^) —^ + {{2I3 + I)x~a)- ^ ' 



dx"^ dx 

+ f /^^ + e„ + + + ^ - - + ^) : -(-^ ' ^ - ] Pi^--^ ix)^0. 

\ 1 + x^ J 

(4.12) 

Making use of the freedom in a and (3, the coefficient in front of 1/(1 + a;^) may nullify, 

a + a^ -b^ + ^ + (3- (3^ = 0, (4.13) 

(1 

-b-2ab+--al3 = 0. (4.14) 

Then, Eq. (|4.12p reduces to the Romanovski equation (|2.18p . Finally, the identification of the 
constants in Eqs. (j4.12p and (|2.18p leads to a condition that defines the energy spectrum of the 
hyperbolic Scarf potential as 

+ en = ~n{2p + n) . (4.15) 
Resolving the three equations (j4.13p . (|4.14p . and (|4.15p for a, /? and e„ results in 

P — —a , a — 2b , 

e„ - -{a-n)\ (4.16) 

This expression for the energy coincides with Eq. (|4.7p . as it should be. In this way it is proved 
that the D'-^'^^x) functions that enter the solution of the Schrodinger equation are equal to the 
Romanovski polynomials. Therefore, the D functions are polynomials. As a result, the wave 
functions in x space take the form 

g„(a;) = (l + a;2)-te-'"^""'(^)£)(,-'^'2b)(^)^ dx = ^/TTx^dz . (4.17) 
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The weight function from which the D polynomials are obtained via the Rodrigues formula is 

^(a+i.-^26)^^^ = (1 +x2)-'^-^e-2'"'^"^'^ (4.18) 
The wave function is now equivalently rewritten to 

* / rfsinh-i(x) 

V dx 

and is of the type in Eq. (|3.13p . As a consequence, the orthogonaUty integral between 
the wave functions will recover the orthogonality between the polynomials as shown 
in Eq. (|4.24p below. In order to relate the a and (3 parameters to those of the Romanovski 
polynomials one can compare the coefficients in front of the first derivatives in the respective 
Eas. ((4T2)) . and ((2l8)) . 

2{-p+l)x + q^{2l3 + l)x~a, (4.20) 

giving 

l3 = -a = -p+-, -a = q = -2b. (4.21) 
In this way, the polynomials that enter the solution of the Schrodinger equation will be 

^(;3=-a.„=2b)(^) ^ ^ (x) . (4.22) 

They are obtained by means of the Rodrigues formula from the weight function w^'^~^'^'~^''\x) as 

In this fashion, the hyperbolic Scarf potential has been solved in terms of the real Romanovski 
polynomials. 

The orthogonality integral of the Schrodinger wave functions gives rise to the following orthog- 
onality integral of the Romanovski polynomials, 

9^{x)g^,{x)dx^ / {l + x')-'^+h~''''^'^''^-'>Ri:^^--''\x)Rl:,^-^--"'\x)dx, (4.24) 

-C30 J —OQ 

which coincides in form with the integral in Eq. (j2.16p and is convergent for n < a. 

In order to relate the result obtained by us to the current literature, it is quite instructive to 
compare Eq. (|2.18p to the Jacobi equation, 

(1 _ " , ^ ' +h-S-ij + S + 2)x)^^V^ -nin + j + 6+ l)P/j'^)(a:) = . (4.25) 

ax-' ax 

Upon complexification of the argument, x ix, the latter equation transforms into 

(1 + " / ^ + z(7 - 5 - *(7 + '5 + 2)x) °^ " ^ + n(n + 7 + <5 + 1)^^^'^) (zx) = . (4.26) 

ax'' ax 

From a formal point of view, Eq. (|4.26p can be made to coincide with Eq. (|2.18p for the following 
parameters: 

7=-p-|, 5-7*. (4.27) 

As long as identical equations have solutions that differ by at most a phase factor, the Romanovski 
polynomials are related to the complex Jacobi polynomials via 

RlP''>\x) = i''p!i'''''^'^''^'^\ix) . (4.28) 
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In this sense one relates in the hterature the Jacobi polynomials of complex arguments and in- 
dices to the solutions of the hyperbolic Scarf potential. However, this relation is in our opinion 
misleading because the finite orthogonality makes the real orthogonal Romanovski polynomials 

to a specie that is fundamentally different from the complex |pi ^ 2 • p+ 2 j,_ 

The orthogonality properties of the complex Jacobi polynomials depend on the interplay between 
the integration contour and the parameter values and need special care [S^. Equation (|4.28p in 
combination with Eqs. (|2.16p . and (|2.17p in fact states that the contour over which Jacobi poly- 
nomials of the type Pn^''^ \ix) are orthogonal is the real axis and not, as one naively would have 
expected, the finite interval 



4.2 Polynomial construction and finite orthogonality. 

The construction of the i?i ^ ^ ' \x) polynomials needed in the exact solutions of Scarf II is 
now straightforward and based upon the Rodrigues representation in Eq. (|2.14p where we plug in 
the weight function from Eq. (|4.18p . In carrying out the differentiations we find the lowest four 
(unnormalized) polynomials as 

(4.29) 
(4.30) 
(4.31) 

+ {-72b + Mab-24:a'^b)x^ + 2{-2 + a){~15+16a-4:a^)x^ . (4.32) 
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Figure 4.3: The wave function -ipilz). 

The finite orthogonality of the above polynomials has been proved among others in Refs. [27l 
I37j . We here rather shall illustrate this property in terms of the polynomial normalization constants 
which becomes especially transparent in the interesting limiting case of the sech^(z) potential (it 
appears in the non-relativistic reduction of the sine-Gordon equation (c.f. [54j )) where one easily 
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Figure 4.4: The wave function ^p2{z). 




Figure 4.5: The wave function t/j^iz). 
finds that the normalization constant, Nn"'''^^'^\ is given by 




a > 1 , 
;-2a)^ a>2, 

- 16a + 15)^ a > 3 etc. (4.33) 



Using symbohc softwares such as Maple and Mathematica is quite useful in verifying the results 
reported here. The latter expressions show that for positive integer values of the a parameter, 
a = n, only the first {n— 1) Romanovski polynomials are orthogonal, as it should be in accord with 
Eq. (|2.17p . and the comment after (|4.7p . The general expressions for the normalization constants 
of any Romanovski polynomial are defined by integrals of the type /^^(l + a;^)"~^'e'*'''^ ^^^''dx 
and are analytic for [n — p) integer or half-integer. 
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4.3 Romanovski polynomials and non— spherical angular func- 
tions. 



The Romanovski differential equation appears in the problem of a particle within a non-central 
scalar potential, a result that can be concluded form Ref. 43 . In denoting such a potential by 
V{r, 9) , one can make for it the specific choice of 

V2{0) 



Vir,0)^Vi{r) + 



V2{e) = -ccot{e). 



(4.34) 



An interesting problem is the electrostatic non-central potential in which case Vi (r) is the Coulomb 
potential. The corresponding Schrodinger equation, 



2n 



1 



,2 1 



]_d_ 2_9 l_ 

^d^^ d^^ r2 sin(6l) 89 ' d9 ' sin^ (9) dip^ 



^■sH9)-^ 



V{r,9) 



is solved in the standard way by separating variables. 



*(r,6',(p) = E^{r,9,(p), 
(4.35) 



*(r,6',(/?) =7^(r)e(6l)e™'^. 
The radial and angular diff'erential equations for TZ{r) and Q{9) are then found as 



d^nir) 2 d7^(r) 



dr^ 



and 



(Pe{9) 



cot{9) 



r dr 
dei9) 



d6'2 ' "^""'"^ de* 
with 1(1 + 1) being the separation constant 



1(1 + I) 



2nV2{9) m' 



sin\9) 



n{r) = 0, 
6(0) = 0, 



(4.36) 

(4.37) 
(4.38) 



4.3.1 Radial equation. 

The radial equation (j4.37p . in the p variable introduced as r = 
states) reads 

2d^ 1(1 + 1) k 1 

J- 5 — n{p) + { 

p ..p 

with 



k EE 



2 I E 



-^^^P and for E<0 (bound 
)7^(p) = 0, (4.39) 

(4.40) 



This differential equation is identical to the radial equation of the Schrodinger equation for the 
hydrogen atom. For this reason, the substitution TZ{p) = p^G{p) seems convenient. It leads to 



P^ + [(2/ + 2)-p]^ + (.-/-l)G(p) = 0, 

which coincides with the associated Laguerre differential equation, 

(PL{x) , ^ , .dL(x) , -r , . 
""^^ + (/^ + 1 - + W = 0' 

upon setting (3 = 21 + 1. For A, the following condition holds valid, 

/ , dLi(x) 1 , ^,d'^a(x) 

^-=-^'-r^^+2("'--^)^ 



(4.41) 



(4.42) 



(4.43) 



which gives = k — I — 1 
energy is obtained as 



n,.. This implies k 



/ + 1, and with k from Eq. (I4.40p . the 



2h? {r 



l + l) 



(4.44) 
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4.3.2 Angular equation. 

Before proceeding further we wish to notice that for V2{0) = 0, and upon changing variables from 
6 to cos (6'), Eq. (|4.38p transforms into the associated Legendre equation and correspondingly 

9(0) ^^(^1^", P,"(cos(0)), 

an observation that will become important below. Following Ref. [43j one begins by substituting 
the polar angle variable by a new variable, z, introduced via 9 = f{z). This leads to the new 
equation 



dz2 



r(f) 



+ /'(z)cot(/(z)) 



d 

d^ 



-^J^V,{fiz)) + l{l + l)~ . f 

sm^(/(z)) 



i^{z) =0, 
(4.45) 

with f'{z) = and 'il){z) defined as ^{z) = 9(/(z)). Next one can require that f'{z) ap- 

proaches zero at z = like sin(2;), meaning, lim^^o /'(-z)/ sin(z) — > 1, and define f{z) via 



r(f) 



-/'(2)cot(/(z)), 



(4.46) 



In considering the left hand side of the latter equation as the logarithmic derivative of f'{z), and 
the right hand side as the derivative of lnsin(/(z)), Eq. (|4.46|) equivalently rewrites to 



leading to 
From this equation one concludes 



(ln/'(z))'=(ln(sin(/(z))))', 
dfiz) 



dz 
dz 



= sin(/(z)). 
dj{z) 



(4.47) 
(4.48) 



sin/(z)' 

The latter expression can be integrated by means of the so called A function and is given by 



(4.49) 



z = In I csc(/(z)) - cot(/(z))| = A (/(z) - |) 

Finally, upon exponentiation one arrives at 

= |csc(/(z))-cot(/(z))| 
1 - cos(/(z)) 



(4.50) 



(4.51) 



sin(/(z)) 

In using sin(/(z)) = 2 sin (^^) cos and 1 - cos(/(z)) = 2 sin^ ' ^Q- (|43T|) simplifies 



to 



tan 



With that the function /(z) is obtained as 

e = f{z) = 2tan"i(e^). 
From this equation one obtains e^, and as 



e" = tan I - I , e 



cotl- 



(4.52) 



(4.53) 



(4.54) 
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and finds cosh(z) as 

1,/ ^ e^ + e" 
cosh(z) = 



tan(^)+cot(^) 



sin(/(z))' 



(4.55) 



Getting back to /(z) = 9 results in 



1 



Correspondingly, 



cosli(z) = . , sinh(z) = — ■\/cosli2(z) — 1 = — cot(6). (4.56) 
sin(t^) 

cos(6') = -tanh(z). (4.57) 

Now, the derivative of f{z) is calculated as f'{z) — sin(/(z)) — secli(z). Upon substituting the 
last relations into Eqs. (|4.34p . and (j4.45p . one arrives at 



dV(2) 



i^{z) = 0. (4.58) 



dz2 

Defining c = the latter equation takes the form: 

, \ ' + \l(l + l)sech2(z) - 5 sech(z)tanh(z) - m^l Mz) ^ 0. (4.59) 
dz^ '■ ■' 

This equation has same form as the radial Schrodinger equation with Scarf II, i.e. 

—J^ + [i-b^ + a(a + l))sech2(z) - 6(2a + l)sech(z)tanh(z) + e] V(^) = , (4.60) 
provided, 

• l{l + 1) plays the role of — (6^ — a{a + 1)), 

• plays the role of — e, 

• c plays the role of —b{2a + 1). 

Recall that the solution to eq. (j4.60p was obtained in subsection 4.1.1 of the present chapter as, 

i>n{z) = iV„(l + sinh2(z))-te-''*^""(^'"'^(^»i?i°+^'"'"(sinh(z)), 
e„ = —(a — n)^ , — oo < x < +oo , (4-61) 

with Nn being a normalization constant. Back to the 9 variable and in making use of the equality 
sinh(z) = — cot(6'), we find 

e(0) = V„(sinh-i(- cot(e))) = 7V„(l + cot2 (0))-te-^ ™*(^»ii'i"+^'"'''\-cot(0)), (4.62) 

showing that the angular part of the exact solution to the non-central potential under consideration 
is defined by the Romanovski polynomials. The two parameters of the Romanovski polynomials 
have to be determined from the system of three equations, 

-b^ +a{a + l) = l{l + l), (4.63) 

-6(2a + l) = c, (4.64) 
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e„ = —(a — n)^ = —rn^, a > n, m > 0, (4.65) 

meaning that the I, m, and c constants can not be independent. There exist various choices for a 
and b. If defined on the basis of the first two equations, one encounters 

1 

+ 2 




= o - '+o +\ l+o +^ ■ (4-66) 



Substitution of a into the third equation imposes a constraint on Z as a function of m, c, and n. A 
second choice for a and b is obtained by expressing a from the third equation in terms of m, and 
n as a = m + n and substituting in the second equation to obtain b as 

b = -— . (4.67) 

Then the first equation imposes the following restriction on I 



X'^^ib'-aia + l)), l^-l + ^^+X. (4.68) 

This / value which is not necessarily integer, can be plugged into Eq. (|4.44p leading to a (discrete) 
spectrum that no longer bears any resemblance to the 0(4) degeneracy. This is the path pursued 
by Ref. [43j. I here instead take a third chance and express a, b, and c as functions of I alone 
according to 

a = b = l{l + l), n = a-m^l{l + l)-m, c = -6(2a+l). (4.69) 

This choice allows to consider integer I values. In making use of Eqs. (|4.2ip . (|4.22p . the polar angle 
part of the wave function in this case becomes 

(4.70) 

The complete angular wave function now can be labeled by / and m (as a tribute to the spherical 
harmonics) and is given by 

Zr(^,V') =V'„=iO+i)-™(sinh-i(-cot(0)))e""'^ = 

(1 + cot2(0))-^e-'('+i)*-"(--*(^))i?J,%'lt?"''^'^'^Vcot(0))e^^^^ . (4.71) 

It reduces to the spherical harmonics Yi"^{9,(p) for a = 6 = 0. In this way, the Romanovski 
polynomials shape the angular part of the wave function in the problem under consideration. In 
the following, we shall refer to Zl"{9,(p) as "non-spherical angular functions". In the appendix 
I present a side by side comparison of | Yi"^{d,ip) \ and | Zl^{9,(p) \. A comment is in order 
on \ZJ"{9,(p)\. In that regard, it is important to become aware of the fact already mentioned 
above that the Scarf II potential possesses su(l,l) as a potential algebra, a result reported by 
Refs. [TTl[60j among others. There, it was pointed out that the respective Hamiltonian, H, equals 
H — —C— J, with C being the su(l, 1) Casimir operator, whose eigenvalues in our convention are 
j{j — 1) with j > versus j{j 4-1) and j < in the convention of [TTl [S^. As a consequence, the 
Scarf II solutions can be viewed as representation spaces of irreducible SU{1, 1) representations. 
Specifically, in the case under consideration the represenation is dicrete, unitary and of infinite 

dimensionality. It is the one denoted by {Dj^ (™ ^ {9, ip)}, with m' = j, j + 1, j + 2, The SU (1, 1) 

labels m', and j are mapped onto ours via 

m' = a+i = Z(Z + l) + i, j = m'-n, m' = j, j + 1, j + 2, .... (4.72) 
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meaning that both j and m' are a half-integer. In terms of these labels the energy rewrites as 
e„ = — (j/ — i)^. The condition a > n translates now a.s j > In result, 0(6') becomes 



6(6*) = i^n=m'-j (Sinh"l(-C0t6l)) = ^{1+ C0t2 6))-'"' + 5e-26 tan-1 (- cot 6) Jli^'^-^b) ^_ ^ 

(4.73) 

Here we kept the parameter 6 general because its value does not affect the SU {1,1} symme- 
try. Within this context, the ji/'m'-j (sinh~"'^(— cot0) |'s can be viewed as absolute values of 
{Dj^ '™ ^{0, (p)} eigenvector components [9] and realized in terms of the Romanovski polynomials. 
The \ZJ^{9, ip)\ functions are then images in polar coordinate space of {D^_ ^ -'C+^^+s) 

J — Tn.-\- 2 

eigenvector components. The representations are infinite because for a fixed j value, to' is bound 
from below to toJj^jjj — j, but it is not bound from above. For example, |Z|(0, (p)| refers to 

D+_, i'^'=^){e,ip), \Z^{0,ip)\ refers to D+_, (" -^)(0,^) etc. 

J~ 2 J~ 2 

4.3.3 Romanovski polynomials and associated Legendre functions. 

Next, it is quite instructive to consider the case of a vanishing ¥2(0), i.e. c — 0, and compare 
Eq. (|4.58p to Eq. (|4.60p for b = 0. In this case, and in accordance with Eq. (|4.2ip 

l = a = p-^, m^ = {l-nf, q=-2b=0, (4.74) 

which allows one to relate n to I and m as m — I ~ n. As long as the two equations are equiv- 
alent, their solutions differ at most by a constant factor. This allows to establish a relationship 
between the associated Legendre functions and the Scarf II wave functions. In taking into account 
Eqs. (j4.11[) .and (|4.17|) together with Eqs. (|4.56p . one finds cot{9) = — sinh(z) which produces 
the following new relationship between the associated Legendre functions and the Romanovski 
polynomials 

F,"(cos(6l)) - (1 + cot2(6l))-5i?|i+J'°^(- cot(6')) , / - m = n = 0, 1, 2, ...I. (4.75) 

In substituting the latter expression into the orthogonality integral between the associated Leg- 
endre functions, 

P,'"(cos(6i))P,r (cos(6i))dcos(6i) = , l^l', (4.76) 

-1 

results in the following integral 

V + cot2(0))-^i?^it^'°\-cot(0))i?^/lt^^'"^-cot(^))dcos(^) =0, l^l'. (4.77) 
When rewritten to conventional notations, the latter expression becomes 

X = sinh(2) , l-n^l' -n' = to > . (4.78) 

This integral describes orthogonality between an infinite set of Romanovski polynomials with 
different polynomial parameters (they would define wave functions of states bound in different 
potentials). This new orthogonality relationship does not contradict the finite orthogonality in 
Eq. (j2.17p which is valid for states belonging to same potential [equal polynomial parameters). 
Rather, for different potentials, Eq. (j2.17p can be fulfilled for an infinite number of states. To see 
this let us consider, for simplicity , n = n' = I ~ m, i.e. I — V . Given p = I + the condition in 
Eq. (|2.17p defines normalizability and takes the form 

2(Z-m) < 2(; + i) - 1 = 2Z, (4.79) 
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which is automatically fulfilled for any m > 0. The presence of the additional factor of (1 + x'^)~^ 
guarantees convergence also for m = 0. Equation (j4.78p reveals that for parameters attached to the 
degree of the polynomial, an infinite number of Romanovski polynomials can appear orthogonal, 
although not precisely with respect to the weight function that defines their Rodrigues representa- 
tion. The study presented here is kindred to Ref. and Eq. (|3.86p from above. Also there, the 
exact solutions of the Schrodinger equation with the trigonometric Rosen-Morse potential (em- 
ployed as quark-di-quark interaction) have been expressed in terms of Romanovski polynomials 
(not identified as such at that time) and also with parameters that depended on the degree of the 
polynomial. Also in this case, the n-dependence of the parameters, and the corresponding varying 
weight function allowed to fulfill Eq. (|2.17p for infinitely many polynomials. 

4.4 Hyperbolic Scarf potential in the Klein-Gordon equa- 
tion. 

The final example to be considered is the case of the relativistic Klein-Gordon equation, 

{P^ - = , P,=id,, (4.80) 

where /i is the mass. One can introduce two different potentials in this equation. The first is a 
vector potential, A^^ introduced via minimal coupling as P^, — gA^, with g being a constant, and 
the second is a scalar potential, S, introduced via m — > fi + S. The vector potential in the so 
called Coulomb gauge satisfies W ■ A — but on many occasions one simplifies the problem in 
choosing A = in which case only the time like component of the vector potential, gAo, denoted 
by V in the following, enters the equation. In effect, the Klein-Gordon equation with vector and 
scalar potential (in units of c = ft = 1) takes the form 

[(^^ - ^M)' +^'- (^(^) + ^)']^(^) = 0- (4.81) 

The latter equation simplifies significantly when S and V are equal. It has been shown in Refs. [611 
l62] that in this case the solution of Eq. (|4.8ip can be found from those of an associated Schrodinger 
equation. Indeed, for time-independent potentials, the total wave function can be written as 
^(r,*) = e-'^*ip{r), with E being the relativistic energy. This substitution results in: 

[V2 + (Vir) - - {S{r) + Ai)']V(r-) = 0. (4.82) 

From now onward I shall focus on the special case of equal scalar and vector potentials, i.e. 
V{r) = S{r). In the following we change variable to V{r) . In this case, the Klein-Gordon 

equation then rewrites to 

[V^+(i^-^)'-(^+.)W) = 0. (4.83) 

In the following v{r) is taken as the central hyperbolic Scarf potential. Separating variables in 
polar coordinates, ip{r) = R{r)H{9)K{ip), leads to 

K"{ip) + m^K{ip) = 0, (4.84) 

H"{e) + cot{e)H'(e) - [m^ csc^{e) - s{s + l)]H{e) = O, (4.85) 

and 

(r^R'ir))' - [s{s + 1) + {E + ^i)r^v{r) - [E"^ - fi^y]R{r) = 0, (4.86) 

where s and m are the separation constants. The angular equation for the aziniuthal coordinate, tp, 
has solutions satisfying periodical conditions: K{(p) — -^e""'^, m = 0,±1,±2, ... The equation 
for the polar angle, 9, can be transformed by means of a; = cos{6) and becomes 

r,^(Pf(x) df(x) r ITl^ . „, ^ / N 
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It can be identified with the associated Legendre diS^erential equation whose solutions are f{cos{9)) = 
P™(cos(0)). The radial equation can be solved for s = 0. For R{r) = D{r)/r it takes the form 

- ((-B + /i) ^2 + {E + ^l) {B^ -A^- A)sech^ (r) 

+ {E + fi)B{2A + l)sech(r)tanh(r) - {E + n){E - fi)j D{r) = 0. 

(4.88) 

Upon naming 



(E + ^-l)i-A^ + E - ^j.) as e, 



{E + n)A^ as a? 



• {E + - A^ - A) as h"^ - - a, 

• {E + fi)B{2A + 1) as b{2a + 1) , 

and in changing variables to x = sinh(r), D{r) — > f{x), amounts to 

,1 + + + f -'-' + °(° + ') + dei±i), + ^ = 0. ,4.8.,) 

da;^ da; \ i + a;^ 1 + / 

The latter equation is equal to the polynomial form of the Id-Schrodinger equation for the Scarf 
II in Eq. (|4.10p whose solutions have been explicitly constructed in Eq. (|4.19p above. Matching 
parameters leads to 



a = 26, /3 = -a, e„ = -(a-n)2. (4.90) 

The energies are then found as 



1 _ ^2 + 2An - 2A^fi + y/A^ + AA^n - + 4A(A + 2n)n + 4^2 

^" ~ 2(1 + ^2) ' 

2 ^2 + 2An - 2A^^ - ^A^ + 4A3n - 4^2 + 4A(A + 2n)/i + 4^2 

= 2{iTA^) ■ (^-^^^ 

The two values for the energies correspond to particles, and antiparticles, as expected from the 
relativistic Klein-Gordon equation. 
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Figure 4.6: The non-central potential V{r,9), for c = —5, here displayed in its intersection with 
the a; = plane, i.e. for r = y'^ + z^, and 9 = tan^^ ^. The polar angle part of its exact 
solutions is expressed in terms of the Romanovski polynomials. 
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Chapter 5 

Conclusions and outlooks. 



In this thesis I presented 

• the classification of the (orthogonal) polynomial solutions to the generalized hypergeometric 
equation in the respective schemes of Koepf-Masjed-Jamei [5S], Nikiforov-Uvarov [S], and 
Bochner [47] . 

• the explicit construction of the Hermite, Laguerre, Jacobi, Bessel and Romanovski polyno- 
mials, 

• the solutions of the Schrodinger equations with the respective one- and three dimensional 
oscillator, the Coulomb- and the hyperbolic Rosen-Morse potentials, and the ^ e"^ barrier 
in terms of one of the above polynomials. 

As new results I report 

• the exact wave functions of the bound states within the hyperbolic Scarf potential in terms 
of the Romanovski polynomials, 

• the finite orthogonality of the Romanovski polynomials in terms of their normalization con- 
stants, a property that allowed to map the finite number of bound states within Scarf II 
onto a finite set of polynomials, 

• the Romanovski polynomials as main designers of non-spherical angular functions of a new 
type, which we identified with components of the eigenvectors of the infinite discrete unitary 
SU(1,1) representation, {D+_ (™ '='('+1)+^) (0, ^)}, 

• a non-linear relationship between Romanovski polynomials with parameters attached to 
the degree of the polynomial and the associated Legendre functions which lead to a new 
orthogonality integral for an infinite series of such Romanovski polynomials, 

• the solution of the Klein-Gordon equation with equal scalar and vector potentials, taken as 
the hyperbolic Scarf potential. 

I conclude that the Romanovski polynomials represent the most adequate degrees of freedom in 
the mathematics of the hyperbolic Scarf potential. 

Further conclusions are: 

• The orthogonality integral of Schrodinger wave function of the form given in Eq. (|3.13p and 
free parameters always recovers the orthogonality of the involved polynomials. 

• Equation (j3.13p did not hold valid on several occasions in which the parameters happened 
to depend on the degree of the polynomials, in which case the orthogonality of the wave 
functions failed in recovering the orthogonality of the polynomials with the free parameters. 
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• As a rule, polynomials with parameters running with the degree of the polynomial, appear 
orthogonal with respect to a weight function which is altered in comparison with the one 
that enters the Rodrigues formula. 

• The case of the hydrogen atom was special in so far as there the anomalous orthogonality 
integral between Laguerre polynomials was observed for polynomial parameters which did 
not depend on the degree of the polynomial , however the variable x did. 

In future research one can 

• exploit the relationship between the Romanovski polynomials and the associated Legendre 
functions to write down various new recurrence relations for the former, 

• employ their weight function as an extension of the student's t distribution and test it in 
statistical problems of estimating standard deviations from data, 

• study symmetry relationships between the non-spherical angular functions which lead to 
equality between such functions of different parameters as visible by inspection from the 
Appendix. 

The Romanovski polynomials are interesting mathematical entities in their own and future research 
is expected to shed more light on their properties and physics applications. 
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Chapter 6 

Appendix. 

In this Appendix I present a graphical side by side comparison of the absolute values of the 
spherical harmonics | Y™{9,(f) \ and the | Z™{9,(p) \ functions. 




Figure 6.1: | Yo^{e,^) \ vs | Zl^{0,^) \ 
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Figure 6.2: | Y°{9,(p) \ vs | Z°{e,ip) \ 
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Figure 6.6: \ Y^^^^O , ^) \ vs \ Z^\0 , ^) \ 
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Figure 6.9: \Yf'ie,^)\vs\Zi\9,^)\ 
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Figure 6.10: | ¥^"^^{0,^) \ vs | Z^^ie,^) \ 
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Figure 6.11: | Y^{9,ip) \ vs | Z2{0,(p) \ 
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Figure 6.17: | Y^\e,^) \ vs | Z^\e,^) \ 
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Figure 6.22: | Zi{9,(p) |, and | ^{9,if) \ as examples for non-spherical angular functions 
don't have a spherical counterpart. 
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